Comparative analysis of the differences between using LiDAR and contour-based DEMs for hydrological modeling of runoff generating debris flows in the Dolomites

Present work aims to explore the differences in hydrological modeling when using digital elevation models (DEMs) generated by points from LiDAR surveys and those digitized on the contour lines of the regional technical map (RTM) and their relevance for the simulation of debris flow triggering. Hydrological models for mountainous areas are usually based on digital elevation models (DEMs). DEMs are used to determine the flow path from each pixel, by which the basin is discretized, to the outlet. Hydrological simulations of runoff that triggered debris flows occurred in two rocky headwater basins of Dolomites, Fiames Dimai (area approximately 0.03 km2) and Cancia (area approximately 0.7 km2) are carried out using a DEM-based model designed for simulating runoff that descends from headwater areas. For each basin, the runoff is simulated using DEMs that are generated using points from LiDAR, and those digitized on the contour lines of the regional technical map, respectively. The results show that the peak discharge values corresponding to the simulations carried out using the LiDAR-based DEMs are higher than those corresponding to the simulations carried out using the RTM-based DEMs. Larger differences are observed for the Dimai basin, where the area corresponding to the RTM-based DEM is markedly smaller than the area corresponding to LiDAR-based DEM, whereas for the Cancia basin, the two areas are similar. Both the differences in the peak discharge and the basin area are due to the poor accuracy of the contour-based DEM (i.e., elevation accuracy), that is, a poor representation of the morphological features that leads to errors on the watershed divide and simplifications of the flow paths from each cell to the outlet. This result is highly relevant for estimating the triggering conditions of runoff generated debris flows. An incorrect simulated value of peak discharge can lead to errors both in planning countermeasures against debris flows and in predicting their occurrence

Evaluating the Differences of Gridding Techniques for Digital Elevation Models Generation and Their Influence on the Modeling of Stony Debris Flows Routing: A Case Study From Rovina di Cancia Basin (North-Eastern Italian Alps)

Debris \ufb02ows are among the most hazardous phenomena in mountain areas. To cope with debris \ufb02ow hazard, it is common to delineate the risk-prone areas through routing models. The most important input to debris \ufb02ow routing models are the topographic data, usually in the form of Digital Elevation Models (DEMs). The quality of DEMs depends on the accuracy, density, and spatial distribution of the sampled points; on the characteristics of the surface; and on the applied gridding methodology. Therefore, the choice of the interpolation method affects the realistic representation of the channel and fan morphology, and thus potentially the debris \ufb02ow routing modeling outcomes. In this paper, we initially investigate the performance of common interpolation methods (i.e., linear triangulation, natural neighbor, nearest neighbor, Inverse Distance to a Power, ANUDEM, Radial Basis Functions, and ordinary kriging) in building DEMs with the complex topography of a debris \ufb02ow channel located in the Venetian Dolomites (North-eastern Italian Alps), by using small footprint full- waveform Light Detection And Ranging (LiDAR) data. The investigation is carried out through a combination of statistical analysis of vertical accuracy, algorithm robustness, and spatial clustering of vertical errors, and multi-criteria shape reliability assessment. After that, we examine the in\ufb02uence of the tested interpolation algorithms on the performance of a Geographic Information System (GIS)-based cell model for simulating stony debris \ufb02ows routing. In detail, we investigate both the correlation between the DEMs heights uncertainty resulting from the gridding procedure and that on the corresponding simulated erosion/deposition depths, both the effect of interpolation algorithms on simulated areas, erosion and deposition volumes, solid-liquid discharges, and channel morphology after the event. The comparison among the tested interpolation methods highlights that the ANUDEM and ordinary kriging algorithms are not suitable for building DEMs with complex topography. Conversely, the linear triangulation, the natural neighbor algorithm, and the thin-plate spline plus tension and completely regularized spline functions ensure the best trade-off among accuracy and shape reliability. Anyway, the evaluation of the effects of gridding techniques on debris \ufb02ow routing modeling reveals that the choice of the interpolation algorithm does not signi\ufb01cantly affect the model outcomes

Modeling Shallow Water Flows on General Terrains

A formulation of the shallow water equations adapted to general complex terrains is proposed. Its derivation starts from the observation that the typical approach of depth integrating the Navier-Stokes equations along the direction of gravity forces is not exact in the general case of a tilted curved bottom. We claim that an integration path that better adapts to the shallow water hypotheses follows the "cross-flow" surface, i.e., a surface that is normal to the velocity field at any point of the domain. Because of the implicitness of this definition, we approximate this "cross-flow" path by performing depth integration along a local direction normal to the bottom surface, and propose a rigorous derivation of this approximation and its numerical solution as an essential step for the future development of the full "cross-flow" integration procedure. We start by defining a local coordinate system, anchored on the bottom surface to derive a covariant form of the Navier-Stokes equations. Depth integration along the local normals yields a covariant version of the shallow water equations, which is characterized by flux functions and source terms that vary in space because of the surface metric coefficients and related derivatives. The proposed model is discretized with a first order FORCE-type Godunov Finite Volume scheme that allows implementation of spatially variable fluxes. We investigate the validity of our SW model and the effects of the bottom geometry by means of three synthetic test cases that exhibit non negligible slopes and surface curvatures. The results show the importance of taking into consideration bottom geometry even for relatively mild and slowly varying curvatures

The debris flow occurred at ru secco creek, venetian dolomites, on 4 august 2015: Analysis of the phenomenon, its characteristics and reproduction by models

On 4 August 2015, a very high intensity storm, 31.5 mm in 20 min (94.5 mm/h), hit the massif of Mount Antelao on the Venetian Dolomites triggering three stony debris \ufb02ows characterized by high magnitude. Two of them occurred in the historical sites of Rovina di Cancia and Rudan Creek and were stopped by the retaining works upstream the inhabited areas, while the third routed along the Ru Secco Creek and progressively reached the resort area and the village of San Vito di Cadore, causing fatalities and damages. The main triggering factor of the Ru Secco debris \ufb02ow was a large rock collapse on the northern cliffs of Mount Antelao occurred the previous autumn. The fallen debris material deposited on the Vallon d\u2019Antrimoia inclined plateau at the base of the collapsed cliffs and, below it, on the Ru Salvela Creek, covering it from the head to the con\ufb02uence with the Ru Secco Creek. The abundant runoff, caused by the high intensity rainfall on 4 August 2015, entrained about 52,500 m3 of the debris material laying on the Vallon d\u2019Antrimoia forming a debris \ufb02ow surge that hit and eroded the debris deposit covering the downstream Ru Salvela Creek, increasing its volume, about 110,000 m3 of mobilized sediments. This debris \ufb02ow routed downstream the con\ufb02uence, \ufb02ooding the parking of a resort area where three people died, and reached the village downstream damaging some buildings. A geomorphological analysis was initially carried out after surveying the whole basin. All liquid and solid-liquid contributions to the phenomenon were recognized together with the areas subjected to erosion and deposition. The elaboration of pre and post-event topographical surveys provided the map of deposition-erosion depths. Using the rainfall estimated by weather radar and corrected by the nearest rain gauge, about 0.8 km far, we estimated runoff by using a rainfall-runoff model designed for the headwater rocky basins of Dolomites. A triggering model provided the debris \ufb02ow hydrographs in the initiation areas, after using the simulated runoff. The initial solid-liquid surge hydrographs were, then, routed downstream by means of a cell model. The comparison between the simulated and estimated deposition-erosion pattern resulted satisfactory. The results of the simulation captured, in fact, the main features of the occurred phenomenon

Coarse-grained debris flow dynamics on erodible beds.

A systematic set of flume experiments is used to investigate the features of velocity profiles within the body of coarse-grained debris flows and the dependence of the transport sediment concentration on the relevant parameters (runoff discharge, bed slope, grain size and form). The flows are generated in a 10 m long laboratory flume, initially filled with a layer consisting of loose debris. After saturation, a prescribed water discharge is suddenly supplied over the granular the runoff triggers a debris flow wave that reaches nearly steady conditions. Three types of material have been used in the tests: gravel with mean grain size of 3 and 5 mm, and 3 mm glass spheres. Measured parameters included: triggering water discharge, volumetric sediment discharges, sediment concentration, flow depth and velocity profiles. The dynamic similarity with full-sized debris flows is discussed on the basis of the relevant dimensionless parameters. Concentration data highlight the dependence on the slope angle and the importance of the quasi-static friction angle. The effects of flow rheology on the shape of velocity profiles are analyzed with attention to the role of different stress generating mechanisms. A remarkable collapse of the dimensionless profiles is obtained by scaling the debris flow velocity with the runoff velocity, and a power law characterization is proposed following a heuristic approach. The shape of the profiles suggests a smooth transition between the different rheological regimes (collisional and frictional) that establish in the upper and lower regions of the flow, and is compatible with the presence of multiple length scales

Atmospheric circulation patterns, cloud-to-ground lightning, and locally intense convective rainfall associated with debris flow initiation in the Dolomite Alps of northeastern Italy

The Dolomite Alps of northeastern Italy experience debris flows with great frequency during the summer months. An ample supply of unconsolidated material on steep slopes and a summer season climate regime characterized by recurrent thunderstorms combine to produce an abundance of these destructive hydro-geologic events. In the past, debris flow events have been studied primarily in the context of their geologic and geomorphic characteristics. The atmospheric contribution to these mass-wasting events has been limited to recording rainfall and developing intensity thresholds for debris mobilization. This study aims to expand the examination of atmospheric processes that preceded both locally intense convective rainfall (LICR) and debris flows in the Dolomite region. 500 hPa pressure level plots of geopotential heights were constructed for a period of 3 days prior to debris flow events to gain insight into the synoptic-scale processes which provide an environment conducive to LICR in the Dolomites. Cloud-to-ground (CG) lightning flash data recorded at the meso-scale were incorporated to assess the convective environment proximal to debris flow source regions. Twelve events were analyzed and from this analysis three common synoptic-scale circulation patterns were identified. Evaluation of CG flashes at smaller spatial and temporal scales illustrated that convective processes vary in their production of CF flashes (total number) and the spatial distribution of flashes can also be quite different between events over longer periods. During the 60 min interval immediately preceding debris flow a majority of cases exhibited spatial and temporal colocation of LICR and CG flashes. Also a number of CG flash parameters were found to be significantly correlated to rainfall intensity prior to debris flow initiation

Design of two retention basins along the torrent Liera on the Gares Valley (Dolomites, North East Italy) after the storm Vaia

A storm called “Vaia” affected the North East Italy and South West of Austria at the end of October 2018. On the Gares Valley (Canale d’Agordo, North East Italian Alps), the abundant runoff triggered several in-channel debris flows that transported hundreds of thousands of cubic meters of sediments on the valley bottom. Some control works for retaining the sediment volume have then been built in the most threatened sites. Herein we show the design of two retention basins that should protect both the Liera torrent and the main road parallel to it. We computed the debris-flow volume to be retained both using an empirical law and simulated solid-liquid hydrographs corresponding to a return period of 300 years. The estimate of the debris-flow volume in the empirical law depends on the basin area upstream of the deposition zone, whereas, in the solid-liquid hydrographs, it depends on the area of the basin closed at the initiation area, as well as on the design rainfall